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C^ . We present a macroscopic model of mixed multi-lane freeway traffic that can 
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I. INTRODUCTION 

Simulating traffic flow is not only of practical importance for developing traffic optimiza- 
tion measures [jl|. It is also interesting because of the various self-organization phenomena 
resulting from non-linear vehicle interactions 0-§]. This includes the formation of traffic 
jams or stop-and-go waves [§|-0], and of synchronized congested traffic ||8|-[T0|. Many of these 



phenomena can be simulated by one-dimensional models like cellular automata ITT.T^ and 
macroscopic traffic models P,p!^. It is, however, questionable if traffic dynamics at on- and 
off-ramps or intersections can be correctly described without consideration of lane changes 



1^-|19|- The same holds close to changes in the number of lanes, accidents, or construction 
sites. It is, therefore, desireable to have multi-lane models that explicitly take into account 
overtaking maneuvers and lane changes. In contrast to existing cellular automata models 



for multi-lane traffic P0|-p3|, we will focus on macroscopic models, here, since they allow 
analytical investigations and a simple calibration to empirical data. 

Apart from lane-changing maneuvers, traffic dynamics is considerably influenced by the 
composition of traffic into various types of vehicles with different desired velocities and 
acceleration capabilities. This can even cause new kinds of phase transitions in mixed 



traffic, e.g. to a coherent, solid-like state of motion [^. Thus, it is also favourable to have 
a macroscopic model that allows to distinguish several vehicle types. 

We will obtain such a model from a generalized version of a gas-kinetic traffic model, from 
which we managed to derive a one-dimensional model that is consistent with all presently 
known properties of traffic flow [0, including synchronized congested flow ||lOl. Although 



lane-changes and several vehicle types make the model quite complex, it is still possible to 
evaluate the Boltzmann-like interaction terms. Moreover, whereas previous approaches have 
neglected correlations between successive vehicles, we will also show how these can be taken 
into account. 

Our calculations take care of the fact that vehicles do not interact locally, but with 
the next vehicle in front, so that the interaction point is advanced by about the safe vehicle 
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distance. Nevertheless, we were able to evaluate the interaction integral without the gradient 
expansion used in previous publications p^j25[| . This resulted in characteristic non-local 



and anisotropic interation terms, which have very favourable numerical properties |2^. Like 
diffusion or viscosity terms, they cause a smoothing of shock fronts. However, they do 
not change the character of the set of partial differential equations from a hyperbolic to a 
parabolic one. Finally, the non-local terms allow to write the macroscopic traffic equations in 
flux form with source/sink terms, so that numerical standard procedures can be used for their 
robust and efficient numerical integration [^ . As a consequence, it is now possible to solve 



the multi-lane traffic equations without elimination of the dynamic velocity equations, i.e. 
without averaging over self-organized velocity variations like stop-and-go traffic. This had 
to be done in previous approaches because of numerical stability [^, with the consequence 
that the investigation of unstable multi-lane traffic was not possible. Now, we are able to 
study how the vehicle dynamics on one lane influences the others. 

Our paper is organized as follows: In Section II we will introduce the kinetic model, 
which allows the derivation of the macroscopic multi-lane model of heterogeneous traffic 
presented in Section III. The decisive steps of this derivation are given in the Appendix. 
Section IV will discuss the calibration of the model to real traffic data for the case of one 
effective vehicle type, and it will display various simulation results for difficult test scenarios 
like lane closures. A comparison with results of an effective single-lane model will be the 
topic of Section V, while Section VI summarizes the results of this paper. 

II. THE KINETIC MULTI-LANE MODEL 



Prigogine and coauthors [|28| were the first who proposed a kinetic equation for the 
phase-space density of vehicles on a highway. Later, Paveri-Fontana ^^ made important 
improvements, formulating an equation for the "extended" phase-space density g{x, v, vq, t), 
where x is the location, v the actual velocity, and vq the desired velocity of a vehicle at time 
t. Nagatani and Helbing have suggested extensions to multi-lane traffic [^,^. However, 
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the validity of all theses equations was restricted to the regime of small densities. 

Here, we present the general kinetic model for multi-lane traffic composed of different 
types of spatially extended vehicles, and derive the corresponding macroscopic model from 
it. The main features of the model are 

• the correct description of dense traffic by taking into account the finite space require- 
ments of vehicles, 

• the "non-locality of the interactions" , which means that drivers look ahead and adopt 
their behavior to the traffic situation at some position in front of them, 

• the replacement of the "vehicular chaos" assumption used by most previous kinetic 
traffic models p4| , p5| , p8| , |29| , p| (with a few exceptions [^,0) by much less restrictive 



assumptions, which account for possible velocity correlations of interacting vehicles. 

For simplicity, we assume that all vehicles of the same type a G {1, . . . , A} have the same 
desired velocity Vq° that may depend on the lane i G {!,..., A^}. Then, we can represent 
the "extended" phase-space density gi{x,v,VQ,t) of all vehicles in each lane through a sum 
of phase-space densities /"(x, v, t) of vehicles of different types a: 

g^ix, V, vo, t) = E /r(x, V, t)5{v, - V,^) . (1) 

a 

The phase-space densities obey the following continuity equation 0,^: 

The terms on the left-hand side represent the continuous change of the phase-space density 
due to the movement of a point in phase space. The third term represents the tendency of 
drivers to accelerate to their desired velocity Vq" with a relaxation time r^. The terms on 
the right-hand side represent (quasi-)discontinuous changes of the phase-space densities due 
to lane-changing maneuvers and braking interactions between vehicles. 
We consider three different types of lane-changing behaviour: 
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• lane-changing maneuvers due to interactions between vehicles, which correspond to 
overtaking maneuvers or lane changes in order to avoid intensive interactions with 
slower vehicles in a lane, 

• "spontaneous" lane-changing maneuvers, which reflect effects of traffic regulations and 
the internal tendency of drivers to use a particular lane, 

• "mandatory" lane changes due to lane mergings, on- or off-ramps, accidents etc. 
Accordingly, we have: 

Following Prigogine, we speak of "interaction" , when a faster vehicle approaches a slower 
one, which forces the faster vehicle to change lane or brake. For simplicity, we will assume 
that vehicles do not change their velocities during lane-changing maneuvers and that braking 
vehicles decelerate exactly to the velocity of the vehicle in front. 

The key quantities, which determine the RHS of the kinetic equations, are the "interac- 
tion rates" , or the effective number of interactions between vehicles per unit of time. Let us 
denote by X"^(a;, f ,t) the interaction rate of vehicles of type a and velocity v at place x in 
lane i with slower vehicles of type b and velocity w < v in front; by J'°'^{x, v, t) the number 
of interactions in lane i between vehicles of type a with velocity w > v and vehicles of type b 
in front with velocity v. The interaction rates Tf'{x, v, t) with 6 G {!,..., A} contribute to 
the decrease of the phase-space density f°'{x,v,t) due to lane changing or braking, whereas 
J'°'^{x,v,t) contribute to the increase of the phase-space density due to braking of vehicles 
of type a to velocity v. They are given by the following formulas: 

Xf = X^ I dwiv- w)ft{x, v,x + St, w, t) , (4) 

Jt = X^ I dwiw- v)ftix, w,x + st, V, t) . (5) 

w>v 

That is, the interaction rates are proportional 
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• to the "effective cross section" Xi = Xi{Pi): reffecting the increased number of interac- 
tions in dense traffic due to the vehicular space requirements (see below), 

• to the relative velocities |f — w| of interacting vehicles, and 

• to the pair distribution function f^^{x, v, x+s^, w, t), describing the phase space density 
of having in lane i two vehicles of types a and b at places x and xf = (x + s^) with 
velocities v and w, respectively. Note that the classical "vehicular chaos" assumption 
would correspond to the factorization /f ''(x, v,x + s°, w, t) = /f (x, v, t)f^{x + s", w, t) 
of the pair distribution function into the phase space densities of single vehicles. 

With s" > we assume an anticipative driver behavior, i.e. a reaction to the traffic situation 
at the advanced "interaction point" xf = (x + s^). A reasonable formula is s^ = 7j"(l/p™'^^ + 
T^Vi), where pf"^^ is the maximum vehicle density in lane i, T" is the safe time headway of 
vehicles of type a, and T^Vi is the safety distance at an average velocity of Vi in lane i. For 
j"- = 1, the vehicles react to the traffic situation at the safe vehicle distance, for jf > 1, 
they look further ahead. The anticipation factor 7" is typically between 1.0 and 3.0 [T^l. 



Now, let us denote by p" the probability that a vehicle of type a in lane i can change lane 
without any delay. Evidently, this probability equals the sum of probabilities of changing 
to the right lane {i — 1) or to the left lane (i + 1): p'^ = p"j„i + vli+v We assume these 
probabilities to be functions of the macroscopic variables such as vehicle densities. 

In terms of the interaction rates, we can now specify the rate of braking interactions by 



\ dt 
and the lane-changing term due to interactions by 



1-Pt)j:[^r\x,v,t)-Ifix,v,t)], (6) 



(^) = E \l^-i,^tiix, V, t) - p^.iJf (x, V, t) 

+P?+i,X^i(a;, V, t) - Pl,^,lf{x, V, t)] . (7) 

The rates of spontaneous lane-changes are simply proportional to the phase-space den- 
sities of vehicles. The proportionality factors are determined by the transition rates 1/T°'_^ ■ 
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of this process, where T°_^- are the characteristic times between lane changes. We will 
assume that they are functions of the macroscopic variables like the densities. Thus, the 
spontaneous lane-changing terms have the form: 

\ ^Wlc -^j-l,j -M^i-1 

-'-i+l,i -'-iyi+l 

An appropriate specification of the transition rates is the following: 

This form is rather simple and, at the same time, it is in accordance with empirical data 
1 3^ . The results displayed in this paper are for /3i = and /?2 = 8. 

The form of the mandatory lane-changing term depends on many different factors in- 
cluding visibility or the location of traffic signs indicating the end of a lane. An explicit 
example will be given in Section IV. 

III. THE MACROSCOPIC MULTI-LANE MODEL 

One advantage of the kinetic equation is, that it allows the systematic derivation of 
equations for macroscopic variables. The macroscopic variables of interest are the densities 
p^{x,t), average velocities V"j"(x, t), and velocity variances O^i^x^t) of vehicles of type a in 
lane i at place x and time t. They can be defined as moments of the phase-space densities: 



pt{x,t) = J dvfnx,v,t), (10) 

Vtix, t) = plix, ty' J dv vftix, V, t) , (11) 

9i{x,t) = pt{x,tr' fdv{v-vn'fnx,v,t). (12) 



One can obtain the macroscopic equations from the kinetic equation, using an iterative 
procedure, called Chapman- Enskog expansion PJ53|] . The first step of this procedure gives 



the so-called Euler-like traffic equations, which are in good agreement with empirical findings 
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|10| , p!3[| . To obtain these, one assumes the velocity distribution in every point x to be in a 
local equilibrium. Although the exact equilibrium solution of the kinetic traffic equations is 
not known, the Gaussian distribution of velocities is a good approximation ||3^. Next, recall 
that we need to specify the form of the pair distribution function of interacting vehicles. 



Previous models [^,^,^,^,0 assumed the pair distribution function to be a product 
of Gaussian one-vehicle distribution functions, which corresponds to the "vehicular chaos" 
assumption that interacting vehicles are uncorrelated. Here, we only presuppose that the 
pair distribution function can be approximated by a general bivariate Gaussian distribution 
function: 



/f (x, v,x + si w, t) = p^x, t)p\{x + si t)^^e-^«(^'-) , (13) 

where B{v, w) is a general positive definite quadratic form and detB the determinant of the 
corresponding symmetrical matrix. Thus, we take into account possible correlations between 
the velocities of interacting vehicles. One can express the coefficients of B{v, w) in terms of 
the moments of the distribution, namely the variances 6^ and the correlation coefficient kf'^: 

The associated determinant equals 



-1 



detB={etef[l-{kff]y . (15) 

A prime indicates that the respective quantity is evaluated at the advanced interaction point 
xf = (x + s°) rather than the actual position x. 

In order to evaluate the macroscopic equations, one needs to multiply the kinetic equa- 
tion by v'', integrate over the velocity v, and close the system of equations by a suitable 
approximation for the moments of higher order [^,^,^ . We will derive the Euler-like equa- 
tions for the vehicle densities and average velocities, and close the system by approximations 
for the variances and correlation coefficients. 

The details of the calculation are given in the Appendix, here we present only the final 
results. 
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The density equations are 



d d 



A 



je{J-l,i+l}6=l 



ab\ 



je{j-i,i+i} V i.* 



'T'a 



and the equations for the traffic flows read: 

= §-s^s. - yn - (1 - p") E Bfibvt) 



h=\ 



+ E E h-.q'(^v;" 

iG{i-l,J+l}6=l 

pv^t p\yr 



]^,cf{svr 



+ E 



j 3 



'T'a 



(16) 



(17) 



One can obtain the corresponding velocity equations by inserting (p!6D into (p^T]). However, 
we will use the flux equations (|T^ instead, because they are more suitable for numerical 
integration methods. 

The functions Af', Bf', and Cf* in the above equations are denoted as Boltzmann factors, 
since they originate from the Boltzmann-like interaction integrals @) and (|^). Thus they 
describe the influence of interactions on traffic dynamics. The Boltzmann factors A'}'^ deter- 
mine the lane-changing flows due to interactions in the density equations, Bf'' the braking 
term, and Cf'' the lane-changing terms due to interactions in the flow equations. The exact 
form of these terms is as follows (for brevity we omit indices of lanes and vehicle types) : 



A{6V) = x{p)ppWS[N{6V) + 6VE{6V)] , 
B{6V) = x{p)ppS \6VN{6V) + (1 + 6V^)E{6V) 
C{6V) = x{p)ppS 



''-Msv)U'-'^ 



V 



6V E{6V) 



,VS ' ' \ S ,/S 

Here, N{z) is the standard Gaussian distribution, and E{z) denotes the error function: 



(18) 
(19) 
(20) 
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N{z) = ^ , (21) 

E{z) = ^ f dye-"^. (22) 

— oo 

The values SVf''' are dimensionless velocity differences between interacting vehicles, which 
are defined by 

^yab ^ y^ ^ ^ (23) 



where 



Sf = 91 - 2kf^9'i9\ + 91 . (24) 

The Boltzmann factors are negligible for negative velocity differences, and they grow rapidly 
with increasing positive differences. The reason for this is rather intuitive: Faster vehicles 
in front do not influence vehicles at the given place x, while slower vehicles force them to 
brake or change lane. 

The dimensionless values (^) have the meaning of "effective" velocity differences. Ac- 
cording to their definition, they increase with 

• the increase of the absolute velocity difference, 

• the decrease of the variance, and 

• the increase of the correlation coefficient. 

The last two properties require some explanation. The increase of the variance as well as the 
decrease of the correlation coefficient both lead to an increase of the factor 5*. This results in 
the decrease of the effective velocity difference and, consequently, of the dimensionless parts 
of the Boltzmann terms (standing in square brackets). Nevertheless, the Boltzmann factors 
themselves increase (see Fig. |1|, which shows the dependence of Boltzmann factor B on the 
absolute velocity difference for different values of S) . This effect has a clear interpretation: 
The dimensionless parts of the Boltzmann factors describe the influence of the difference in 
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the average velocities at location x and the advanced interaction point x' on the number of 
interactions. An increase of S reduces the effect of this difference in comparison with the high 
variation of individual velocities in the vehicle flow. Thus, the increase of 5* diminishes the 
value of the effective velocity difference 5V , although it enlarges the interaction rates. One 
could also say: An increase of 5* causes higher interaction rates, but also a wider transition 
region between the limiting cases {Vl" — V^'' ) ^ and {V^ — V^* ) ^ 0. 
For the variances ^", we use constitutive relations of the form 

et = at{p.)Vt\ (25) 

according to which the variance of vehicle velocities is a certain proportion a" of the squared 
average velocity, which depends on the total vehicle density 

P^ = T.Pt (26) 

a=l 

in the respective lane i. This is well justifled by empirical flndings. The appropriate expres- 
sion of the functions (y.1{pi) is given below. 

The correlation coefficients could be approximated as functions of the densities at the 
points X and xf = (x + s"), as well as the distance sf: kf' = k^^{p1,p\ ,sf). However, an 
empirical determination of this function is very difficult, as it requires a thorough analysis of 
a huge amount of single-vehicle data. Therefore, we will apply the common approximation 
k"^^ K. for the time being. 

Another important function to be estimated is the "effective cross section" Xi = Xi{pi)- 
This value reflects the increase of the effective number of interactions in dense traffic. In a 



previous publication on the single-lane variant of the above model |]T3| , it was shown that 
the following expression for the effective cross section is consistent with the limiting cases 
at high and low vehicle densities and well justifled by the resulting properties of the model: 

[1 - p{p)]x{p) = -^r--.T, — 7 — V2 ■ (27) 

Note that, without further assumptions, x{p) is determined in the single-lane model only 
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together with the overtaking probabihty p{p). In our multi-lane model, we suggest the 
following decomposition of the above expression: 

/ N _ 1 I K)i^i Pi ,^o\ 

a. . _ exp (-pg.p./p— ) 

'^'^'^'^~ UPi) • ^ ^ 

Here, variables without an index for the specific vehicle type a represent weighted averages 
of the variables belonging to the different vehicle types in lane i, for example, 

T, = Y^f^Tr (30) 

a=l pi 



IV. CALIBRATION OF THE MULTI-LANE MODEL AND SIMULATION 

RESULTS 

Next, we present the results of calibration and simulation for a special case of the general 
model discussed above. We consider the two-lane variant of the model for a single vehicle 
type, where i = 1 represents the right ("slow") lane and i = 2 the left ("fast") lane. The 
calibration was done on the basis of empirical data for the Dutch two-lane highway A9. 

The variance prefactor or "structure factor" ai{pi) = 6i/Vi^ can be estimated on the 
basis of direct observation. The empirical data show a "step-like" but smooth dependence 
of this prefactor on density, with an increase at about 40 vehicles per kilometer. It can be 
well fitted by the following function (see Fig. H): 

1 -1 



ttj(Pj) = Oioi + Aoi 



II P^- Pc 

1 + exp 



(31) 



Spi , 

This step-like form plays an important role, as it determines the specific shape of the 
equilibrium velocity-density relation V^{pi) (see Fig. ^ and the fundamental diagram 
Qti.Pi) = Pi^iipi)- Other parameters that influence the fundamental diagram are the safe 
time headway Tj, the desired velocity Voj, and the maximum density p™^^. The desired 
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velocity fits the inaxiinuin velocity in free traffic, whereas the safe time headway and the 
maximum density determine the slope of the fundamental diagram at high densities and 
its intersection point with the density axis. The acceleration relaxation times Tj and the 
anticipation factors 7^ do not influence the fundamental diagram. Instead, they allow to fit 
the stability behavior and dynamics of traffic flow [|13|. 

The parameters that affect lane-changing processes are the coefficients gij, f3i, and /?2 for 
the spontaneous lane-changing rates, and the coefficients poi for the overtaking probabilities. 
The main sources for estimating these parameters are empirical data on lane occupancies 
and lane-changing rates as a function of density. One can easily obtain the data on lane 
occupancies from the usual measurements made by induction loops. In contrast, the direct 
measurement of lane-changing events is much more difficult, which results in a lack of reliable 



data in the literature ||3^. Luckily, this kind of data is only necessary for the estimation 
of the order of magnitude of the model coefficients, while the ratios of the coefficients for 
different lanes, which essentially deflne the multi-lane dynamics, can be well estimated by 
the available lane occupancy data. 

Figures ^ through ^ show the corresponding flts. Following Ref. 1^^, we assume that 



the maximum lane changing rate of about 500 to 550 events per hour, kilometer and lane 
is achieved at densities of about 20 to 25 vehicles per kilometer. The higher occupancy of 
the left lane at middle and high densities yields higher estimated values of the coefficients 
Qij, poi for the left lane relative to those for the right lane. The primary use of the right 
lane at small densities (see Fig. ^ reflects the European traffic regulations. One can take 
this into account by a "European-rules" correction prefactor of the spontaneous coefficients: 
gi2 -^ gi2gEu{pi),g2i -^ g2i/gEn{pi), where < gEnipo) < l is a smooth step-like function 
similar to (|3T|). 

The calibration results show that the spontaneous lane-changing terms influence mainly 
the low-density regime, while the lane changes due to interactions, which are negligible at 
small densities, determine the difference in lane occupancy at high densities. For the typical 
flt of the lane occupancy curves (see Fig. ||) it turned out that the overall contributions 
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of spontaneous and interactive lane-changing terms to the total lane-changing rate were 
approximately the same (Fig. H). Note that, at low densities, lane changes correspond 
mainly to interactive lane changes from the right to the left lane and to spontaneous lane 
changes from the left to the right lane, which is plausible for European traffic. 

The results of the parameter estimation are summarized in Table |. 

One important property of the above model is the ability to describe the development 
of different congested traffic states |38|j3^ . Figure |^ shows the development of stop-and-go 
traffic, which arises from a small density perturbation in the right lane. Due to lane changes, 
the perturbation spreads to the other lane, and the traffic dynamics on the neighboring lanes 
becomes synchronized |]8|j9|j40[| . In particular, this holds for the propagation of large density 
clusters. Nevertheless, the traffic flow in the left lane behaves more unstable in the range 
of moderate densities. This fact is in agreement with observations and can be theoretically 
explained by the different velocity-density relations (decreasing more rapidly for the left 
lane, see Fig. ^. 

Next, we present multi-lane simulations of the interesting case of a bottleneck, corre- 
sponding to an on-ramp, a lane closure, or an accident. Assume, for example, that the right 
lane ends at certain place Xend- We expect that the resulting traffic situation will depend on 
the volume of the incoming flow. 

In order to model the behavior of traffic close to a bottleneck, we must specify the 
mandatory lane changes. In the framework of the above multi-lane model, this can be 
done by introduction of additional lane-changing terms describing a sufficient increase of 
lane changes to the left lane, whereas lane changes to the right lane will be surpressed by 
setting the corresponding coefficients close to the bottleneck to zero. However, the following 
difficulty arises: While the density on the right lane decreases to zero at the bottleneck, the 
velocity (which depends on the density and velocity on the neighboring lane) can stay large 
up to the very end of the lane. This causes numerical problems in keeping the density and 
flow positive everywere. To avoid this, we apply the following calculation procedure to the 
last section of the right lane in front of the bottleneck, where we have assumed that this 



Shvetsov/Helbing: Macroscopic Dynamics of Multi-Lane Traffic 15 

section is of length Lq = 500 m throughout this paper. 

Close to the bottleneck, all drivers in the right lane (producing a traffic flow of volume 
piVi) must merge into the adjacent lane. This implies that the drivers in the right lane 
will adopt their velocity Vi{x,t) to the velocity V2{x,t) in the left lane. In addition, we will 
assume that the lane-changing rate grows inversely proportional to the remaining distance 
L{x) = (xend — x), in order to guarantee that all vehicles have changed lane at the place 
Xend where the right lane ends. Hence, for x G [xcnd — Lq, Xend], our model for the right lane 
reads: 

dpi djpiVi) ^ _ piK .^21 

dt dx L{x) ' ^ ^ 

Vi{x,t) = V2ix,t). (33) 

For the left lane, we have: 

dp2 d{p2V2) ^ piVi_ . . 

dt dx L{x) ■ ^ '' 

In order to describe a smooth transition from the "normal" to this "adaptive" behaviour in 
the merging zone, we evaluate the RHS of the equation according to 

RHS = [1 - A;(x)]RHSnorm + A;(x)RHSadapt , (35) 

where k{x) is a smooth step-like function similar to (|31|) with A;(xend~-^o) ~ and k{xend) ~ 
1. 



The results of our simulations are presented in Figures |g through |T0|. The traffic dynamics 
is essentially characterized by the volume of approaching traffic and the capacity of the 
bottleneck (which is given by the outflow from traffic jams in the left lane). At low upstream 
densities (see Fig. §), the capacity of the left lane is sufficient to transport the vehicle flow 
from both lanes. In contrast, we have an immediate formation of congested traffic upstream 
of the bottleneck, if the total traffic volume in both lanes exceeds the capacity of the left 
lane (cf. Fig. |]). Surprisingly, for a certain range of moderate densities, the resulting traffic 
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situation turns out to depend on the initial condition. While a perfectly homogenous flow 
will lead to an increased but free traffic flow downstream of the bottleneck, as for small traffic 
volumes, a small perturbation can trigger the breakdown of traffic flow, although the left 
lane could carry the total vehicle flow in both lanes (cf. Fig. |10|). This can happen, when the 
traffic flow downstream of the bottleneck is unstable. As long as the perturbation is small, 
it moves downstream. However, when its amplitude becomes larger, it eventually changes 
its propagation speed and flnally travels upstream, until it reaches the bottleneck. Then, 
traffic breaks down, and a steadily growing region of congested traffic develops upstream of 
the bottleneck, whereas traffic downstream of the bottleneck flows freely (Fig. |l^). A similar 
phase transition from free to congested traffic is known to occur close to on-ramps ||T0|,|9[|. 



V. COMPARISON WITH THE EFFECTIVE SINGLE-LANE MODEL 

The traffic situations discussed above can be also simulated with an "effective" single- 
lane model that implicitly averages over the dynamics of all lanes. The corresponding model 
was proposed in [|13)|13 ^^"^ basically corresponds to our multi-lane model, applied to one 
lane only, so that the lane changing terms drop out. For our simulations, we use the following 
"effective" model parameters: Vq = llOkm/h, pmax = 150 vehicles/km, r = 35 s, T = 1.6 s, 
7 = 1.2, ao = 0.007, Aa = 0.031, pc = 0.28p"^^^, and 5p = 0.025p'"^^. 



In Figure |Tl], we compare the average of the densities in the left and the right lane 
according to the multi-lane model (see Fig. |^) with the effective single-lane model. It turns 
out that, despite of the sensitive dynamics in the unstable traffic regime, both models 
produce similar spatio-temporal traffic patterns, but there are some differences in detail. 
This shows that the effective single-lane model gives already a reasonable representation of 
the traffic dynamics, although it produces small deviations from the djTiamics predicted by 
the multi-lane model. 

Let us make a similar investigation for the example of a bottleneck. In this case, we 
can treat the merging lanes in the effective single-lane model by a reduction of the effective 
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lane number I{x) from 2 to 1 within the merging section, i.e. for x G [xend ~ -^Oi^^end]- 
For example, we may use the linear relation I{x) = [1 + L{x)/Lo]. The conservation of 
the number of vehicles implies the virtual ramp flow u = —{pV/I) dl/dx, which gives the 
following continuity equation for the effective vehicle density per available lane: 

dt dx I{x) dx ' 



The result of the corresponding simulation is presented in Figure |12| in comparison with 
the plot of the average density per lane obtained with the multi-lane model. The pictures 
show a good correspondence between the multi-lane and the effective single-lane models. 
However, there are slight differences in the form and propagation velocity of the upstream 
front of the congested traffic region. These originate from the fact that the average dynamics 
of two nonlinearly behaving systems with different parameters cannot simply be represented 
by one system of the same type with suitably averaged parameters, as it can be done for 
linear systems. 

VI. SUMMARY AND CONCLUSIONS 

We have proposed a gas-kinetic traffic model for heterogeneous multi-lane traffic and 
systematically derived the corresponding macroscopic traffic model. Thus, effects of different 
vehicle types and lane changes are explicitly taken into account. Whereas previous multi- 
lane models have usually assumed spontaneous lane changes only, we managed to calculate 
the lane-changing rates due to vehicle interactions and found that these are of the same 
order of magnitude. Note that both, spontaneous and interactive lane changes are necessary 
to describe the empirically observed density-dependence of the total lane changing rates, 
the lane occupancies, and the density difference among lanes correctly. 

Moreover, the multi-lane traffic model formulated above treats vehicular space require- 
ments and high vehicle densities in the right way, and we have even discussed possible effects 
of velocity correlations of interacting cars, which basically reduce the interaction rates. The 
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corresponding computer simulations are robust also in the unstable traffic regime, so that 
we did not need to eliminate the dynamic velocity equation, as was done in a previous study. 

We have successfully calibrated our model to empirical traffic data. The resulting model 
is in good agreement with the observed variance-density relations, the velocity-density re- 
lations, and the occupancies of the different lanes, as well as with the density-dependence 
of the lane-changing rates and the density difference among lanes. We were able to show 
the synchronization effect among lanes due to lane changes and could describe the traffic 
dynamics at bottlenecks. A comparison of the average dynamics in the different lanes with 
corresponding simulation results of an effective single-lane model showed a qualitative, but 
not fully quantitative agreement. 

Our present investigations focus on the empirical evaluation of velocity correlations be- 
tween interacting vehicles and on the calibration of the model to a mixture of vehicle types 
like cars and trucks, both of which are difficult tasks. We expect that this will allow us to 



describe the effects of heterogeneous traffic which were found in microscopic models |23,41 



and other approaches P2 
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APPENDIX: 

In order to evaluate the macroscopic equations, one needs to calculate the ffist two 
moments of the kinetic equation in velocity space. This procedure was described in details 
in P, p^ , |35| . The new contribution of this paper consists in evaluating the general Boltzmann 
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factors ([ISD through (^0]) including vehicular space requirements and possible velocity corre- 
lations of successive vehicles. The Boltzniann factors are defined by the following integrals 
of interaction rates: 



Af = Jdvlf{x,v,t), (Al) 

Bf = Jdv vlf{x, v,t)- J dv vjf{x, V, t) , (A2) 

Cf = f dvvlf{x,v,t). (A3) 



Here, we present the evaluation of the integrals of Xf* only, but the integration of J'"''' is 
analogous. For brevity, we will omit the indices of lane and vehicle types in the following. 
This means that the pair distribution function f{x,v,x + s,w,t) actually denotes the pair 
distribution function f^^{x,v,x + s^,w,t) for particular types of vehicles a at point x and 
6 at (x + s") in a particular lane i. Consequently, V, 9 stand for V^{x), 0^{x), and V , & 
stand for Vi{x + s°), B\{x + sf). Also, we drop the multipliers x{p)pp' ■ According to the 
definition of interaction rates (|^), one must evaluate the integrals 

Bk = dv / dw v^{v — w)f{x,v,x + s,w,t) (A4) 



v>w 



for k G {0, 1}. 

From the mathematical perspective, the task is to integrate a bivariate Gaussian distri- 
bution, multiplied by polynomials in v and w, over the half-plane v > w. This can be done 
by linear transformations in the {v, w)-plane in three steps: 

• Transform the bivariate distribution to the canonical rotation-symmetric form, 

• rotate the plane to make the boundary of the integration area parallel to one of the 
coordinate axes, and 

• separate variables. 

Then, the integration over one of the axes becomes trivial, and the integration over the other 
axis gives combinations of terms which contain Gaussian and error functions. 
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The bivariate Gaussian distribution / is defined by ( [I3D and (|T4D . One can transform tfie 
quadratic function (|1^) to the simplest symmetric form vl + vol by hnear transformation: 



(A (v\ 



yojj 



\^'J 



Ci 



f.,\ 



yw,j 



( 1 



Ci 



1 



V V^ 



COS(y9 

siny? 



= sin 09 ' 
A ^ 



cosip 



(A5) 



Here, A± are (positive) eigenvalues of the quadratic form (|T^: 

1 



A- 



2e^'ri - m 



' ±^9"^- 2(1 - 2k'^)ee' + ^'2 , 



(A6) 



and the angle v? is defined by 

1 



tany? = ^- (e' -e- Je^ - 2(1 - 2k^)e9' + ^'2) 



(A7) 



After this transformation, we obtain 



1 



f{x,v,x + s,w,t) dvdw = — e aK+^'i) dy^ dyj^ 

for the pair distribution function, and the boundary of the integration area becomes 

cos 09 — sin 09 cos o? + sin 09 , 

V — w = 7== f 1 7=^ Wi + V — V = . 



Ou 



(A8) 



(A9) 



Next, we apply additional rotation, which does not change the symmetric form of the distri- 
bution, and make the boundary of integration area parallel to one of the axes, say y. This 
is done by 



f.^ 



v^v 



Co 



( \ 

X 



\yj 



Co 



( I ■ l\ 

COS ip , — sm ip 



sin Ip , cos Ip 



(AlO) 



where 



sin Lp — cos Lp cos 09 + sin 09 

cos?/^ = , sm?/; = 



v^sa: 



VS\- 



(All) 



s 



[sin(/9 — cosv9)2 (cos(y9 + sin(y9)' 



A, 



A_ 



e-2kvee' + e'. 



(A12) 
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Finally, we consider the composition of two linear transformations: C = C2 ° Ci = {cij). 
This transformation brings the integral ( [A4| ) into the form with separating variables: 



Bk= I dx'-^ f dy%L^{5V -x){y + cux + ci2yf- (A13) 

J V 27r J V 27r 

x<5V -00 

The integration over y becomes trivial, now, as it corresponds to evaluating the moments 
of a normal distribution. The evaluation for /c G {0, 1} involves only the first two moments, 
which equal to 1 and 0, hence the integration over y in those cases results just in the 
elimination of the integral over y and y-containing terms from the expression above. 
The coefficient cn reads 

cos 09 cos ^ sin 09 sin 'j/' k\fdd' — Q , , . .^ 



Hence, for /c G {0, 1} we obtain 

k 



Bk= J dx^=S{5V-x)i^ + x\ . (A15) 

Note that the integration of JT", which is neccesary to obtain the Boltzmann factor 5, 
leads to the same expression with the last factor under the integral being replaced by 



-"^^'^'-xt 



s 
The remaining task is to evaluate the "incomplete moments" of the normal distribution, 

which can be expressed through the normal distributions lSS[x) and error functions E[x) 

(see notations (|2T| ) and (|22|) ). Applying the formulas 

I dxN{x) =E{a), (A16) 



x<a 



dx xN{x) = -N{a), (A17) 



x<a 



dx x^N{x) = -aN{a) + E{a) (A18) 



x<a 



to the integral above, one obtains the desired expressions ( P^ through (^) for the Boltz- 
mann factors. 
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TABLES 



Parameter 




Notation 


Right lane 


Left lane 


Desired Velocity 




Vo 


105 km/h 


123 km/h 


Maximum Density 




max 


150 vehicles/km 


150 vehicles/km 


Relaxation Time 




r 


35 s 


35 s 


Safe Time Headway 




T 


1.7 s 


1.2 s 


Anticipation Factor 




7 


1.2 


1.2 


Coefficients for Variance Approximation 


ao 


0.007 


0.0065 






Aa 


0.03 


0.036 






Pc 


0.275 p'^'^^ 


0.305 p'^'''' 






5p 


0.03 p'"^^ 


0.025 z?"^'^^ 


Coefficient for Overtaking 


Probability 


Po 


17.0 


12.5 


Coefficient for Spontaneous Lane-Changing 


9i,3-i 


75 


28 



TABLE L The estimated parameter values for the two-lane, single vehicle-class model, cali- 
brated to traffic data from the Dutch motorway A9. 
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FIGURES 
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FIG. 1. Boltzmann factor i? as a function of the absolute velocity difference for different 
values of the variance factor S. 



0.08 

0.07 

g 0.06 

^ 0.05 

cu 

^ 0.04 

o 

.2 0.03 

1- 

> 0.02 

0.01 





Structure Factor 



*++^ 






Left Lane * 


^!;r:<:C^* ***** 





10 20 30 40 50 60 70 
Density (vehicles/km) 



80 



FIG. 2. Fit of the density-depenent variance prefactor or "structure factor" ai{pi) specified in 
Eq. (pTI) ( — ) to the empirical data of the relative velocity variance Oi/V^ {+, x). The corresponding 
parameter values are listed in Table |. Note that the deviation of the fit curve for the right lane from 
the empirical data at small densities is probably a consequence of assuming one vehicle type only 
instead of heterogeneous traffic. However, this deviation is of minor importance for the dynamic 
properties of the model, since it is limited to the stable density regime of free traffic flow. 
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Velocity-Density Relations 
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Density (vehicles/km) 

FIG. 3. Equilibrium velocity-density relations of the multi-lane model ( — ) for the parameter 
values specified in Table | in comparison with empirical data from the Dutch highway A9 (-I-, 
x). The symbols represent the (vertical) averages of one- minute data that were evaluated for 14 
successive days. 
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FIG. 4. Theoretical ( — ) and empirical (-I-) difference of the densities in the left and the right 
lane as a function of the lane-averaged density. 
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FIG. 5. Fit of the lane occupancy, i.e. of the percentage of the vehicle density in a lane 
compared to the total density in all lanes. Symbols correspond to empirical data, lines to the 
results of our multi-lane model. The preference for the right lane at small densities comes from 
the European traffic regulations. 
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FIG. 6. Equilibrium lane-changing rates according to our multi-lane model as a function of 
the average density. Note that, at small densities, we have more interactive lane changes from the 
right to the left lane ("overtaking maneuvers") and more spontaneous lane changes from the left 
to the right lane, as expected for Europe. 
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(a) Density in tlie Left Lane (veliicles/km) 
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(c) Lane Changes from Left to Right (events/h/km) 
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FIG. 7. (a) + (b): Simulation of stop-and-go traffic on a circular road in the regime of unstable 
traffic flow, arising from an initial density perturbation in the right lane, which eventually spreads 
to the left lane. (c) + (d): The lane-changing rates have temporary peaks at the locations, where 
the traffic situation in the neighboring lanes evolves differently. This tends to reduce the differences 
among lanes, so that similar spatio-temporal traffic patterns form in both lanes. Consequently, we 
have a synchronization of lanes at medium and high vehicle densities. 
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FIG. 8. Simulation of a closure of the right lane for x > Xcnd = 6 km at time t = Omin. At 
a density of pinjt = 12.8 vehicles per kilometer and lane or lower (here: pinit = 12.6 vehicles/km), 
the capacity of the left lane is large enough to transport the vehicle flow from both lanes, resulting 
in a higher vehicle density in the left lane downstream of the bottleneck, whereas the right lane is 
empty behind the lane closure. 
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FIG. 9. (a)+(b): Closure of the right lane as in Figure ^, but for an average initial density of 
Pinit = 25 vehicles per kilometer and lane. Behind the bottleneck, a growing region of congested 
traffic forms immediately, since the capacity of the left lane is exceeded by the traffic volume in both 
lanes, (c) Note that the effective capacity of the left lane (i.e. of the bottleneck) is considerably 
less than the maximum flow and the flow in the left lane upstream of the jam. This is, because the 
outflow from congested traffic is a self-organized quantity [13|, which is of the order of 2000 vehicles 
per hour, here. The flow per lane in the jammed region is half of this characteristic outflow, (d) 
The step-like structure of congested traffic corresponding to regions of two different densities in (a) 
and (b) is related to a deceleration in two steps (rough braking and fine braking) , when approaching 
a traffic jam from free traffic. This behavior has been also observed in a microscopic traffic model 
|43|. According to an explanation by Ansgar Hennecke, it relates to the pronounced hump of the 
fundamental diagram in the density region between 20 and 50 vehicles per kilometer. Where the 
traffic flow is stable (at densities around 60 vehicles per kilometer and higher), the flow-density 
relation tends to stay close to the fundamental diagram (i.e. the equilihrium flow-density relation). 
In contrast, the dynamic flow-density relation is a self-organized relation in the density regime 
of unstable traffic flow (at densities between about 25 and 55 vehicles per kilometer), connecting 
the stable flow with the self-organized outflow from traffic jams. Hence, we will usually have 
different slopes in the resulting dynamic flow-density relation, corresponding to the propagation 
of congested regimes of different densities with different speeds. This behavior disappears for 
smoother fundamental diagrams, in which the congested part decreases more or less linearly. 
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FIG. 10. Simulation of a closure of the right lane as in Figured with the same initial density 
of Pinit = 12.6 vehicles per kilometer and lane, but with a small perturbation of the traffic flow in 
the left lane. Altough the total flow of 2590 vehicles per hour in both lanes is below the maximum 
possible flow in the left lane of 2630 vehicles per hour, traffic flow eventually breaks down. In other 
words: If all vehicles would use only the left lane, there would be no traffic congestion upstream 
of Xend' The breakdown of traffic is initialized by a perturbation of traffic flow that eventually 
gives rise to a growing region of congested traffic llO] , from which the self-organized outflow is only 



about 2000 vehicles per hour |13]. 
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(b) Effective Single-Lane Density (vehicles/km/lane) 
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FIG. 11. Comparison of (a) the average density according to the multi-lane model and (b) 
the density resulting from an effective single-lane model for the formation of stop-and-go traffic 
presented in Figure |^. 
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FIG. 12. Comparison of (a) the average density according to the multi-lane model and (b) 
the density resulting from an effective single-lane model for the case of a lane closure displayed in 
Figure |9|. 



